The critical amplitude ratio of the susceptibility in the random-site two-dimensional 

Ising model 
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We present a new way of probing the universality class 
of the site-diluted two-dimensional Ising model. We analyse 
Monte Carlo data for the magnetic susceptibility, introducing 
a new fitting procedure in the critical region applicable even 
for a single sample with quenched disorder. This gives us 
the possibility to fit simultaneously the critical exponent, the 
critical amplitude and the sample dependent pseudo-critical 
temperature. The critical amplitude ratio of the magnetic 
susceptibility is seen to be independent of the concentration 
q of the empty sites for all investigated values of q < 0.25. 
At the same time the average effective exponent 7 e / / is found 
to vary with the concentration q, which may be argued to be 
due to logarithmic corrections to the power law of the pure 
system. This corrections are canceled in the susceptibility 
amplitude ratio as predicted by theory. The central charge of 
the corresponding field theory was computed and compared 
well with the theoretical predictions. 



The effect of impurities on the critical behavior is 
one of the central questions of phase transition phenom- 
ena . According to the Harris criterion, weak concen- 
trations of impurities do not affect the critical exponents 
ii du > 2, where d is the system dimension, and v is the 
exponent of correlation length in the pure system || . In 
the two-dimensional Ising model (2DIM) one has dv = 2. 
In other words, the 2DIM is a marginal case, where loga- 
rithmic corrections may become important in the vicinity 
of the phase transition, while critical exponents are not 
changed, as has been found analytically for the specific 
heat and correlation length Q as well as for the magnetic 
susceptibility ||. Indeed, the widely accepted picture is 
that weak bond dilution changes the critical behavior of 
the correlation length of the pure model £ ~ 1/\t\ by a 
logarithmic term into 
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where r = (T — T c )/T is the reduced temperature, with 
T c being the critical temperature and go is a coefficient 
proportional to the strength of disorder (strictly speak- 
ing, go is the central charge of the N = Gross-Neveu 
model related to the 2DIM, see [Q- ||, and g = for 
the pure 2DIM). Similarly, the critical exponent of the 
magnetization M and the magnetic susceptibility \ ex ~ 



ponents remains the same as for the pure 2DIM, with the 
critical behavior being modified by logarithmic factors. 
This prediction for weak bond disorder was confirmed 
numerically in a number of papers ■ 

In the following, we study the universality class of the 
site-diluted Ising model. The phase diagram of this model 
contains two crucial points, the pure Ising fixed point at 
zero concentration of impurities, q — 0, and the percola- 
tion fixed point at q = q c = 0.407254 §|. The site-diluted 
2DIM was investigated numerically by quite a few au- 
thors pJ- in the last few years. The interpretations 
of the results are rather controversial. Some of authors 
claim |lC],[ll|] that the critical exponents, e.g. that of the 
magnetic susceptibility, 7, vary with the impurity con- 
centration q, but that the ratio, say, 7/1/ (where v is 
the critical exponent of the correlation length) does not 
change. Others |l2|jr^ ] concluded that impurities lead, 
for the susceptibility as well as for other quantities, only 
to logarithmic corrections as for weak bond dilution. In 
all cases, the analyses have been performed either above, 
at or below the critical point, T c . 

The main aim of the present study is to identify the 
universality class by computing the critical amplitude ra- 
tio of the magnetic susceptibility x, thereby comparing 
data both below and above T c . 

It is well known that the universality class is not only 
characterized by its critical exponents but also by the 
critical amplitude ratios [p"4|| . For instance, in zero exter- 
nal magnetic field the critical behavior of the magnetic 
susceptibility x is given by x — r r 7 i n the symmetric 
phase and by x — r" M~ 7 in the ordered phase, with T 
and r' being the critical amplitudes and 7 and 7' are the 
critical exponents. T/T' is the critical amplitude ratio. 

Note that the critical amplitude ratio is very often 
quite sensitive to the universality class. The basic idea is 
that if the values of exponents vary from the pure Ising 
model to the percolation ones as reported in jl0],[ll) , then 
one could expect a variation of the critical ratio T/T' 
as well. So, we may expect that the critical ratio will 
change from the value T/T' = 37.69365 known for the 
pure two-dimensional Ising model fll5| to the percolation 
limit value T/T' 170 |fiq| . Clearly, such a variation of 
about four times is much easier to check than the vari- 
ation of the exponents which is reported to be approxi- 
mately 10 percent. 

To investigate the critical behavior of the magnetic 
susceptibility and the corresponding amplitude ratio, we 
carried out extensive simulations of the 2DIM with site 
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dilution. Numerical calculations of the critical ampli- 
tudes of the susceptibility are known to be difficult, even 
for the pure Ising and Potts models JT7|-f2l| due to finite- 
size effects and analytical corrections to scaling. For the 
diluted case, the problem becomes even more delicate be- 
cause of a lack of quantitative knowledge on the critical 
temperature and the possible crossover associated with 
the randomness induced logarithmic factors. 

The crucial technical idea of the present work is a 
heuristic fitting procedure to determine the critical am- 
plitude ratio. 

The paper is organized as follows. First, we will sum- 
marize the main results. Then, we will describe the fit- 
ting procedure and demonstrate how it works for the pure 
model. Then, specific results will be presented and dis- 
cussed. 

The main results of our study can be summarized as 
follows: 

1. New fitting procedure: A sample dependent pseudo- 
critical temperature T*(q; L, i) (i refers to the sample and 
L is the linear dimension of the square lattice) follows 
from fitting the susceptibility in the high and low tem- 
perature phases close to criticality, to x(j) — P f~ le i'f 
and x(f) ~ r' |f| _7c ^, respectively. 7e//>7e// are aver- 
age effective critical exponents determined in a range of 
temperatures, where finite size effects and corrections to 
scaling may be neglected; f = \T — T* (q; L,i)\/T is the 
sample dependent reduced temperature relative to T*. 
The fitting condition is the equality of the effective crit- 
ical exponents "f e ff and leff, which is applicable if the 
corrections to scaling are not large in the critical region 
as is the case for the 2D Ising model. 

2. Universality of the critical amplitude ratio for the 
magnetic susceptibility: We estimated the ratio T/V by 
fitting our Monte Carlo data for lattices with linear di- 
mension L = 256. The ratio seems to remain constant 
for all concentrations of impurities we considered, i.e. for 
q in the interval [0,0.25]. This behavior may be inter- 
preted as a manifestation of the universality class of the 
site-dilute Ising model being independent of the degree 
of dilution. It should be noted that theory predicts the 
cancellation of the logarithmic corrections ||,[l7 18 1 in the 
ratio of high-temperature and low-temperature suscep- 
tibilities and, thus, predicts the universality of critical 
amplitude ratio, at least for small concentration of im- 
purities where DD theory |i) are applicable. 

3. Weak and strong dilution regions: We found the 
value of dilution q* = 0.1 up to which the predictions 
of the weak-dilution theory j§,|]|8) works well could be 
explained as concentration at which two intrinsic lengths 
coincide. The first is the average distance between im- 
purities H U ex exp(l/g). The second is the percolation 
correlation length £ p oc (q c — q)~ 4 ^ 3 characterizing the 
size of the holes formed on the lattice by the diluted 
sites. For the concentration of diluted sites larger than 



q* the theory of weak dilution cannot be applied because 
the size of the conglomerates of the diluted sites becomes 
larger than the average distance between the diluted sites 
and the effective interaction between impurities should be 
taken into account. 

4. Variation of the average effective exponent Jeff-' 
We find that the average effective critical exponent 7 e / f 
varies with the concentration of impurities, similarly to 
previous observations Jl0|-p^[ . The variation may be at- 
tributed to logarithmic corrections, growing linearly with 
the small concentration of dilution up to the concentra- 
tion q* = 0.1. 

5. Central charge <?o of the N = Gross- Neveu model: 
We extract from numerics its dependence upon the con- 
centration of diluted sites and found the values coin- 
cide well with ones predicted by exact expression for the 
charge in terms of the site-diluted model H for the dilu- 
tion concentration q up to about q* =0.1. 

The model: Each site of a square lattice is either oc- 
cupied by an Ising spin, Si = ±1 or not. The frac- 
tion of empty sites is denoted by q. The positions of 
non-occupied sites were generated using the shift-register 
generator with lags (9689,417) which is known to be 
appropriate for selecting randomly lattice sites p2[ . 
At each concentration of empty sites <?, chosen to be 
q = 0,0.03,0.07,0.1,0.15,0.18,0.2,0.22, and 0.25, the 
number of realizations (or samples) in the simulations 
ranged between ten and twenty five. 

It seems obvious that the singular part of the magnetic 
susceptibility stems from fluctuations of spins belonging 
to the percolating cluster of occupied sites p],p3[|. Those 
spins which are disconnected from the percolation cluster 
do not contribute to the critical behavior. Accordingly, 
we took into account only the fluctuations of the spins in 
the largest cluster, thereby reducing the "noise" in the 
susceptibility due to the small clusters (and reducing the 
computation time). 

The 

simulations: For each concentration of empty sites, we 
computed the magnetic susceptibility \(T,q; L,i) in the 
critical region at about forty temperatures. The temper- 
atures were chosen in the interval t t < Jrl < r Q where r r 
is the rounding temperature r r ~ 1/L |24| . above which 
finite-size effects may be neglected, and r a is the reduced 
temperature, above which corrections to scaling become 
important Q . For the pure Ising model, r a may be esti- 
mated from the exact solution l[15f 



k B T X (t) = r |rP (l + e x r + ...) 



(1) 



with e x = 0.07790315; the corrections to scaling, e x r, 
become important only at rather large reduced tempera- 
tures, say, t > r a w 0.13. Presumably, dilution does not 
affect r a significantly Q. 

In the simulations, we used the one-cluster-flip algo- 
rithm p5[, discarding the first 10 4 clusters for thermal 
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equilibration. Totally, 10 5 clusters were generated for 
each sample i, at given values of q and T. 

The fitting procedure: The fitting procedure is based 
on the assumption that the critical exponent of the mag- 
netic susceptibility 7 takes the same value below and 
above the critical temperature. Thence, one may deter- 
mine the sample dependent pseudo-critical temperature 
T*(q; L, i), as described above, with the average effective 
critical exponents, j e ff and 7^ coinciding in the inter- 
val [r r ,r a ]. In this way, we extracted sample dependent 
critical amplitudes T(i) and T'(i) as well as the effective 
critical exponent 7e//(«)- 

Assessment of the accuracy of our approach: To check 
the accuracy of our approach, we studied also the pure 
Ising model, where the critical amplitudes are known ex- 
actly. 

First, we considered the ratio of xi T ) = (M 2 ) — 
((\M\)) , as computed in the simulations, to the singular 
part of the magnetic susceptibility multiplied by the lead- 
ing correction to scaling (i.e. x{ t \L)/t ~ 7 (1 + e x r), see 
Eq. (1)), with the linear dimension L of the square lat- 
tice ranging from 16 to 256. To reduce finite-size effects, 
we set in the simulations (|M|) = in the symmetric, 
high temperature phase, T c > 0. The results are shown 
in Fig |], where the upper solid line corresponds to the 
exact value of critical amplitude above the critical point, 
r = 0.962582, and the lower solid line corresponds to 
T' = 0.025537. Clearly, the critical amplitudes are ap- 
proximated closely in the range r r < |r| < r a , especially 
for L= 256. Note some deviations of the magnetic sus- 
ceptibility from the asymptotically exact value deeply in 
the ordered phase, which could be eliminated taking into 
account the background term Dq = —0.104133... [ p6[ . 
Actually, we have checked the sensitivity of the fitting 
results to the inclusion of this term. One might estimate 
the accuracy of the numerical value of the critical am- 
plitude ratio r/r' (=37.69..) extracted from the Monte 
Carlo data to be about 5 percent. 

Next, with L — 256 and for the different concentrations 
of dilution q, we consider the ratio of the Monte Carlo 
data for the magnetic susceptibility x/x\ averaged over 
the various samples. The temperature scale has been de- 
termined by using T* (q) . Obviously, as shown in Fig ||, 
the ratio is (nearly) constant in the range [r r , r a ]. In that 
temperature range, one may expect that both finite-size 
effects and corrections to scalings are negligible. How- 
ever, crossover terms in the presumed logarithmic cor- 
rections are expected to depend rather sensitively on the 
degree of dilution, q. However, if they are of similar 
nature above and below the critical point, our approach 
circumvents this difficulty, as seems to be consistent with 
the results depicted in that figure. It is a commonly ac- 
cepted picture Jlj],[l8| that the logarithmic corrections 
should canceled in the ratio of the high-temperature and 
low-temperature susceptibilities. Our results clearly sup- 



port such scenario. 

Specific results and discussion: Our results on the crit- 
ical amplitude ratio and the pseudocritical temperature, 
averaging over the sample-dependent results, for L=256 
and various concentrations of dilution, are presented in 
Table H with 1-cr error bars in parentheses. One may em- 
phasize that the pseudocritical temperature computed by 
us differs from the values of the "true" critical tempera- 
ture reported by other authors Jl(],[ll| (usually based on 
finite-size analyses) by less than CH~r), reassuring us of 
the correctness of our approach. It is remarkable, that 
the ratio of critical amplitudes remains constant within 
the error bars, while the critical amplitudes by them- 
selves grow by almost a factor of three, when varying 
the dilution from q= to 0.25. We limit our simulations 
at concentration q = 0.25 because larger dilutions need 
larger system sizes in order to have a reasonably wide 
critical region. Unfortunately, simulations of larger sizes 
is out of computation power at the moment. The (small) 
deviation of the critical amplitude ratio in the pure case 
from the exact value could be reduced by including in 
the fit the background term Dq. The fit in the reduced 
pseudo-critical temperature window 0.015 < f < 0.2 
gives D = -0.07(2) and ratio T/V = 38.7(4) deviates 
only by two-sigma from the exact value. 

Figure || shows the dilution-dependent critical ampli- 
tudes r and r', given in Table as well as the values 
obtained from the fit to the magnetic susceptibility aver- 
aged first over samples. The nice agreement of the results 
shows that the order of the averaging plays no role. The 
statistical errors are slightly lower in the former case, as 
already mentioned by Wiseman and Domany p7| . 

The effective exponent 7 e / / is also given in the Table |. 
Obviously, starting in the pure case, Jeff first increases 
quite rapidly with dilution, but changes only mildly at 
stronger dilution, q 0.2. The initial variation of 7e//> 
at weak dilution, may be explained quantitatively. Ana- 
lytically, the magnetic susceptibility has been calculated 
P8|,p| to have the form 



X(t) =r|r|" 7 / 4 (l + 0.07790315r)(l-.g ln|r|)' /0 (2) 



7/8 



with the coefficient g given [E9J by 



a = -go = - 



Jr ir (1 + V2/tt) 2 1 - q 



4.843- 



1-9 



(3) 



In Fig ||, the coefficient g, from (^), is plotted together 
with the fit of the magnetic susceptibility data, below 
and above the critical point, to (^|). The results are ex- 
pected to agree for weak dilution; indeed, the agreement 
holds up to g=0.1. Otherwise, pronounced deviations 
set in for stronger dilution. In fact, Fig |] demonstrates 
two facts. On one hand, it is a check of the recent an- 
alytic, supposedly exact result by Plechko on the coef- 
ficient g. On the other hand, it is consistent with (||), 
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which, in turn, readily explains the variation of the av- 
erage effective exponent Jeff, computed in the interval 
[r r , T a ], with dilution (see also Ref. Q). Taking the loga- 
rithmic derivative of (|J) and properly choosing the tem- 
perature interval for the averaging ]30[ | of the In r term, 
one can get a linear dependence of the effective exponent 
jeff oc 7/4 (1 + aq) with a coefficient a of the order of 
unity. This coincides quite well with the values of Jeff 
in the Table for q < q* f» 0.1. 

It seems that this value q* « 0.1 could be explained 
at follows. There are two lenghts for the diluted model 
in addition to the two length scales one has in the Ising 
model: namely the correlation length £ oc 1/r and the 
system size L. The first length ij oc exp(—l/g) is de- 
fined by the value at which the term g In |t| in (§) 
becomes of the order of unity. The next length is the per- 
colation correlation length [^| £ p oc (q c — q)~ 4 / 3 . One 
could check that these two lenghts coincide for q w 0.1. 
Physically this means that for q > q* the disorder is no 
weaker than assumed in the theories mentioned above. 

The critical amplitudes are practically the same for the 
weak dilution q < q* . Also the effective exponents are 
visibly modified in this region by logarithmic corrections. 
They start to grow in the region of the "strong" dilution 
q* < q < q c and their ratio seems to remain unchanged. 
Probably, one could expect to see a crossover regime from 
Ising to percolation universality class in the very vicinity 
of the percolation point. 

Conclusions: Of course, reliable Monte Carlo data on 
even larger system sizes may still be desirable in order to 
check the universality of the ratio of critical amplitudes 
for the magnetic susceptibility, as suggested by our study. 
One should also try to analyse systems with even stronger 
dilution, i.e. closer to the percolation limit p3| . Never- 
theless, our data already provide evidence that the two- 
dimensional Ising model with site-dilution is described 
by the same critical exponent (modulo logarithmic cor- 
rections) and the same critical amplitude ratio of the 
magnetic susceptibility as the pure Ising model, in the 
range of dilutions investigated. The critical amplitude 
ratio r/r' is always quite close to the pure Ising value 
37.69 and far away from the percolation value 160 — 200 
(which is known from simulations [^6| , ^8[ ). The small 
apparent variation of the average effective exponent with 
the degree of dilution may be explained as being due to 
logarithmic corrections. 
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TABLE I. Results of the fit to magnetic susceptibility for 
the lattice size L — 256. 



q 


r' 


r 


Jeff 


T c 


r/r' 


exact 


0.02553 


0.96258 


1.75 


2.26918 


37.685 


0.00 


0.0239(8) 


0.960(5) 


1.757(1) 


2.2684(7) 


40.21(20) 


0.03 


0.0237(2) 


0.946(8) 


1.814(3) 


2.1610(2) 


39.96(19) 


0.07 


0.0239(4) 


0.96(2) 


1.885(6) 


2.0146(3) 


39.98(38) 


0.10 


0.0253(9) 


1.00(3) 


1.93(1) 


1.9032(5) 


39.65(59) 


0.15 


0.0326(1) 


1.27(3) 


1.95(1) 


1.7098(4) 


39.12(55) 


0.18 


0.037(2) 


1.42(6) 


1.99(2) 


1.5905(7) 


38.92(85) 


0.20 


0.043(1) 


1.70(5) 


1.97(1) 


1.5103(4) 


39.98(58) 


0.22 


0.052(3) 


1.99(8) 


1.96(2) 


1.426(1) 


39.1(1.4) 


0.25 


0.066(4) 


2.57(9) 


1.97(2) 


1.298(1) 


40.7(1.9) 
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FIG. 1. The ratios of the computed magnetic susceptibility 
divided by the singular part including the leading corrections 
to scaling (exactly known ([!]) r = |i~|~ 1,75 (l+e x )) for the pure 
Ising model in high (r+) and low (r_) temperature phases for 
system sizes L — 16 (open boxes) , 32 (open circles) , 64 (open 
diamonds), 128 (closed boxes), 256 (closed stars). See text 
for the futher details. 
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FIG. 2. The critical amplitude ratio computed as the 
ratio of susceptibility data in the ordered and the symmetric 
phases x( r )/x(~ T ) f° r various concentrations of dilution, q=0 
(open boxes), 0.03 (open circles), 0.07 (open diamonds), 0.10 
(open stars), 0.15 (closed boxes), 0.18 (closed circles), 0.20 
(closed diamonds) and 0.25 (closed (stars) for the system size 
L = 256. 
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FIG. 3. The critical amplitudes V (closed boxes) and T' 
(closed circles) of the magnetic susceptibility as function of 
the concentration of empty sites q. Open signs denote results 
of the fit to the magnetic susceptibility first averaged over 
samples. Horizontal lines refer to the exact values for the 
pure Ising model. 



g 




FIG. 4. The coefficient g = 4ngo as a function of the im- 
purity concentration. The solid line denotes the analytic re- 
sult of Ref. 8. Circles (stars) correspond to fits of Monte 
Carlo data for the susceptibility in the low (high) tempera- 
ture phase. 



